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| Abstract. Exact solutions of the Dirac-Pauli equation for massive neutrino with 

■ anomalous magnetic moment interacting with dense matter and strong electromagnetic 

ON , field are found. The complete system of neutrino wavcfunctions, which show spin 

rotation properties are obtained and their possible applications are discussed. 

(N ; 

C*** ! 1. Introduction 

o ■ 

The neutrino oscillation phenomenon is well established now pQ. For the detailed 
^ ■ information about neutrino flavour oscillations we refer to the pioneer works [21 [3], 
some resent reviews (H [5] and references therein. From the theory of this phenomenon 
it follows that oscillations are possible only when neutrino possesses non- vanishing mass. 
Consequently, a neutrino has non-trivial electromagnetic properties. In particular, 
a Dirac neutrino has to possess a non-zero magnetic moment [6]. In turn, the 
existence of the magnetic moment is associated with so-called neutrino spin oscillations. 
These oscillations are described as neutrino spin precession or, in other words, as 
transitions between left-handed (active) and right-handed (sterile) neutrino states in 
an electromagnetic field. Phenomenology of this effect was widely discussed in the 
literature [3 El El HQl HH H2l H31 HH UHl HE] . 

However in the most popular papers the theory of neutrino spin oscillation was 
based on solving the Cauchy problem for the Schrodinger type equation with an effective 
Hamiltonian. So the states of neutrino were treated as mixed states. The purpose of 
our work is to find a complete system of solutions of the wave equation which accounts 
for the influence of an electromagnetic field and a dense matter on neutrino dynamics. 
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The existence of such system gives a possibility to describe the states of neutrino with 
rotating spin as pure states and to evaluate the probabilities of various processes with 
neutrino in the framework of the Furry picture. 

2. Model 

When the interaction of a massive neutrino possessing an anomalous magnetic moment 
Hq with background fermions is considered to be coherent, its propagation in matter 
and electromagnetic field is described by the Dirac-Pauli equation with the effective 
potential [TTl [T8] . For a moving and polarized matter such an interaction is described 
by the effective four-potential / M . 
The function 

/ 

is a linear combination of the currents 

3 1 ) = {nfU° f ,n f Uf}, (2) 
and the polarizations 

X>, = + (3) 

of background fermions. In these equations nj and (0 |C/| 2 ^ 1) are > respectively, 
the number density and the mean value of the polarization vector of the background 
fermions / in their center-of-mass system. In this reference frame the mean momentum 
of the fermions / is equal to zero, and = {u°j,Uf} denote the four-velocity of the 
reference frame. Summation is carried out over the background fermions /. We use the 
units h = c = 1. 

The expressions for the coefficients p^' 2 ^ depend on the model chosen for neutrino 
interactions. In the framework of the standard model these coefficients calculated in the 
first order of the perturbation theory have the forms [191 EH EH E2] 

pf = V2G F [l ev + T 3 (/) - 2QW sin 2 # w }, pf = V2G F [l ev + T 3 (/) }. (4) 

Here is an electric charge of the fermion, Tg is a third component of weak isospin, 
G-p and #w are th e Fermi constant and the Weinberg angle respectively. For electron 
neutrino interaction with electrons I eu = 1, positrons I ev = —1, otherwise I eu = 0. 

Let us discuss an equation for neutrino (its mass eigenstate) in matter and 
electromagnetic field. In order to avoid correlations between flavour and spin oscillations 
(the Mikheyev-Smirnov-Wolfenstein effect [TTt [T8]) we shall suppose that the effective 
four-potential is the same for all neutrino flavours and consider the matter, in which a 
fraction of electrons is small, i.e. approximately n e = 0. In what follows, we restrict 
our consideration to the case of a homogeneous medium. Then the explicit form of the 
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generalized Dirac-Pauli equation is uniquely determined by the assumptions similar to 
those adopted in [23] : 

id - \f{l + 7 5 ) - -/ioi^fV - m)v = 0. (5) 



Here F^ u is an electromagnetic field tensor, a^ u = ^(7^7^ — 7^7^). The quantities with 
hats denote scalar products of the Dirac matrices with four-vectors: a = 7 M a M . 

Note, when propagation of neutrino through a dense matter and an electromagnetic 
field is studied, we should keep in mind that the matter and the field are mainly 
situated in the same area. Therefore strengthes of electric and magnetic fields and 
average velocities and polarizations of matter should satisfy the self-consistent system 
of equations including the Maxwell equations, the Lorentz equation 

«J = Z-FSu), (6) 

1 TUf J 

and the classical spin evolution Bargmann-Michel-Telegdi (BMT) equation [21] 



\). (7) 



Here a dot denotes the differentiation with respect to the proper time r. 

As already mentioned, we restrict our consideration to constant velocity and 
polarization of matter and constant homogeneous electromagnetic field in equation (jSJ) • 
However, even such a choice imposes certain limitations on F^ u and f 1 . Indeed, function 
f 1 is a linear combination of currents and polarizations of background fermions. Since 
the velocities, polarizations and number densities of different components of matter are 
independent characteristics in the general case, then from the condition / M = const or 
ffj, — it follows that velocities and polarizations of background components should be 
stationary: ii^ = 0, Ay = 0. Using equations (jSJ), © we can obtain the restriction 

F^f u = 0. (8) 

It should be stressed that this condition is a direct consequence of the fact that the 
average velocity and polarization of particles of matter in an external field should satisfy 
the classical equations of motion. 



3. Calculations 

Let us look for a solution of equation (j3J). Since the functions are constant, the 

canonical momentum operator id^ commutes with the Hamiltonian of this equation. 
However, the commonly adopted choice of eigenvalues of this operator as quantum 
numbers is not satisfactory in our case. Particle kinetic momentum components related 
to its group four- velocity by the relation g M = mw M , q 2 = m 2 are more suitable to 
play this role. In the present paper we discuss the solutions described by the quantum 
numbers which can be interpreted as kinetic momentum components. The explicit form 
of a kinetic momentum operator for a particle with spin is not known beforehand, hence, 
in order to find the appropriate solutions, we have to use the correspondence principle. 
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It was shown in [25] that if effects of neutrino weak interactions are taken into 
account, the Lorentz invariant generalization of the BMT equation for a spin vector 
of a neutrino moving with the four-velocity w M has the form 



where 



(9) 



(10) 



Here and further for the simplification of the formulae we include a value of neutrino 
magnetic moment into an electromagnetic field tensor, i.e. ^oF^ u =>■ F^ v . 

We introduce quasi-classical spin wavefunctions, which can be constructed as 
discussed in [261 EZ] • Suppose the Lorentz equation is solved, and the dependence of 
particle coordinates on the proper time is found. Then the BMT equation transforms to 
an ordinary differential equation, whose resolvent determines a one-parametric subgroup 
of the Lorentz group. The quasi-classical spin wavef unction is represented by a spin- 
tensor, whose evolution is determined by the same one-parametric subgroup. In our 
case this spin-tensor represents the Dirac bispinor. 

Following the technique presented in the works [26, 28j we choose the solution of 
equation © in the form 



ty(x) 



- iF ^U(r(x),r = 0)^ Q (x). 



In this formula U(t, tq) is an operator of evolution for a quasi-classical spin wavefunction, 



-iF(x) ■ 



is a phase factor, and 



-iqx i 



l-Co7 5 £o)(g + m)Vo 



(12) 



is a solution of the Dirac equation for a free particle. Here four-vector Sq determines the 
direction of particle polarization, Qo = ±1 is a sign of spin projection on this direction 
and tpo is a constant bispinor normalized by the condition ty Q (x)^o(x) = m/q^. 
In our case operator U(t, r ) obeys the equation 



U(t,t ) 



-^l\fq ~ qf) + — 2 l 5 H^q ul ^q \ U(r, r ), 
4m m A 1 



(13) 



where H^ v = —\e^ vpX Fp\ is a dual electromagnetic field tensor. It is obvious that the 
solution of equation ( fl3i) can be represented as a matrix exponent 



£/(r,r ) = exp <^ i(r - r ) 



l -l\fq-qf)+ ' - :> 



fH^q vlll q 



(14) 



Am m A 

Substitution of the expression ffl~T]) into equation ([5]) taking into account ([TBI leads 
to the relation 



q + dF(x)-±f + ^ 5 f + 7 5 N 



^-Uq-qf) + — 2 H^q vl ^ q 

Am m A 



l -F^a^ - m e- lF ^U(r(x), r )^ (x) = 0, 
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where N p = <9 M r. Since the commutator [q,U] = 0, and the matrix U(t(x),t q ) is 
nondegenerate, then to hold this relation the following condition is required: 



dF{x) 



^-{fq-qf) + — 2 H^q ul ,q 
4m m z 



2 ^ 



0. 



Since the coefficients at linearly independent elements of the Dirac matrix algebra must 
be equal to zero independently, we receive 

d»F(x) = fP/2, (15) 

and the system of equations for determination of the vector N 11 : 

<p*{m - (Nq)) + {N<p)q» = 0, F^q a = -e^ pX N^ p q x , (16) 

where ^ = f> x /2 + H^q u /m. 

It is easily verified that for the compatibility of this system it is necessary 

q,F^H ua q a + mq,F^f u /2 = 0. (17) 

Since q p takes in general an arbitrary value, in order to satisfy (JT7J) identically it is 
necessary to demand the conditions 

1 



F^H av EE 
F^U = 0. 



0. 



(18) 
(19) 



Hence, the solution of equation (jSJ) can take the form (TTTT) only if the tensor F^ u 
is flat, i.e. its second invariant I 2 = \F^ u H pv is equal to zero, and the vector / M 
is its eigenvector corresponding to the zero eigenvalue. It should be noted that the 
antisymmetric tensor has an eigenvector corresponding to zero eigenvalue if and only if 
this tensor is flat. So conditions (fTBj) and (fl9l are not independent, and (fT8|) follows 
directly from (fl~9l) . Thus we come to the condition above (jSJ), which was obtained only 
on the physical attends. 

Let us introduce an orthogonal basis in the Minkowski space 

m 2 H pv H ua q a - <f TV 



n o = <t l m -> n i 



H^q v 



F^q v 



I In 



where M = q fl H tJlU H up q p , M = q ll F ilv F V pq p . Expanding the four-vector iW with the 
basis and using equations ffl~6]) we get as the result of calculations 



„ A ( PV)W-M-mf l tHi»q t , 

2VM{{^q) 2 - mV) 



1" 



1 + 



m 2 (fip) 



2(((pq) 2 — vrfiip 1 



With the help of the relations 

2(<pq)H^ = m(^f L 



ftp"), Af -Af = 2m 2 h, 



((f V ) 2 - f 2 V 2 )m 2 = (LH^ qi/ ) 2 ) = 8(w) 2 h 



(20) 



(21) 
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where 1% = \F liV F ia/ is the first invariant of the tensor F^ u , expression (1201) can be 
represented as 

^ = _ a M H/V) + f M^L_ + ^ m3 (^) (22) 

y 2((y3g) 2 -mV) J 2(<pq) * 2(ipq)((ipq) 2 - m 2 <p 2 )' K ' 



Let us note that, in order to derive formulae of the form (I2ip . one should take the 
following arguments into consideration. For an arbitrary antisymmetric tensor A pv , its 
dual tensor *A^ V = — ^e flupX A p x, and any four-vectors g p ,h p , (gh) ^ the following 
relation takes place 

A^(gh) = - [g p A vp h p - A pp h p g v ] + * [h»*A V(, g p - *A pp g p h v \ . (23) 

It leads to 

gi *AfAPg»h 2 + (g^ v h v f - h p A»APh»g 2 - (g p A p h»f = 

( 24 ) 

= g*A>fAP v h v {gh) - g p A p A p v h v (gh) = 2(ghfl x . 

Formulae f[2"Tj) follow from equation (|23J) if conditions (IT5|) . ffT9~j) are taken into account. 
Since / M , iV M = const we obtain for the proper time 

r = (Nx), (25) 

and for the phase factor which determines energy shift of neutrino in matter we have 

F(x) = (fx)/2. (26) 

It is not difficult to verify that in view of the relations obtained, the expression for the 
wavefunction takes the form 

= \ E e~ l(Pcx) (l " C7%)(1 - Co7 5 ^o)(g + m)Vo, (27) 



where 



2 

C=±i 



s v = q p (^q)/m- 



= q n + f»/2- (N p ^(tpq) 2 - ¥ 2 m 2 /m = qP(l + C 



2^J ((fq) 2 — m 2 (p' 1 
,1 f „( 1 _ (V(W) 2 ~m 2 ip 2 \ _ C(f<p)m 2 
2 V (cpq) J * 2(<pq)y/(<pq) 2 -m 2 <p* 



(29) 



It is obvious that the system of solutions (|27|) is a complete system of solutions of 
equation (jSJ), which is characterized by the kinetic momentum of the particle q p and 
the quantum number = ±1, which can be interpreted as the neutrino spin projection 
on the direction Sq at the moment r = (Nx) = 0. In the general case this system is 
not stationary. The received solutions are stationary only when Sq = S^ p . In this case 
the wavefunctions are the eigenfunctions of the spin projection operator to the direction 
St p with the eigenvalues ( = ±1 and of the canonical momentum operator id p with the 



Neutrino spin rotation in dense matter and electromagnetic field 7 

eigenvalues P^. The orthonormalized system of the stationary solutions of equation (151) 
can be written down in the following way: 

* ( (x) = e-« p ^^\J\{l - Cl 5 S tp )(q + m)^o, (30) 

where J is the transition Jacobian between the variables and P^\ 



J = det (My) = det 



dPl dPl dq c 



Qqi QqO Qqj 

Using formulae (j2T]) it is possible to rewrite the matrix My in the form 

(/¥>) 
2yJ (<pq) 2 — m 2 ip 2 



Mij = Si, I 1 + C:r^===== I + 



/ nf_\ ( _ m 2 (fip) 2 + f 2 ^q) 2 -m 2 p^ 2 

V 1 ^ (<pq) J Y j 9j Q° J *(<pq)((<pq) 2 - mV 2 ) 3/2 ' 
It can be easily shown that for any vectors g and h the relation 

det (Sij + gihj) = 1 + (gh) 

is correct. Hence 

f ^2a/ (pq) 2 — m 2 (p^j f ^ a/ (pq) 2 — m 2 ip 2 ^ ^ ^ 

The structure of solution ( 1301) directly leads us to the conclusion that when neutrinos 
move through a dense matter and an electromagnetic field which satisfy condition ([19"]) , 
they can behave as free particles, i.e. move with the constant group velocity 
dP? q 

v - = a£ = ? (32) 

and conserve the polarization. 

However in interactions with other particles the channels of reactions which are 
closed for a free neutrino can be opened (see, for example, [2E1 EH]), as a result of 
difference of the dispersion law for the free neutrino P 2 = m 2 and the one for neutrino 
in matter and electromagnetic field: 

P 2 = m 2 - / 2 /4 - 2h - 2CAy / (P$) 2 - <l 2 m 2 , (33) 

where 

pn = pg - ff/2, & = /72 + H^P v /m, 

A = sign (l + c U H ^/ 2m ~ 2I A (34) 
y A/(v?g) 2 -mV J 

In order to get equation (|33|) we used the following relations obtained from formulae 
and 



(Pf) = 2(<pq) 1 + C 



UH^q v /2m - 2l x 



V (v 9 *?) 2 — rn 2 (p 2 
(Pf)(fcp) = 2( W )(/$), $ 2 (^) 2 = ^(^ 
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The appearance of the factor A in equation (1331) is connected with the fact that ( is 
a particle spin projection on the direction determined by a kinetic particle momentum 
instead of the canonical one. 

Note that the solutions which are classified by a kinetic momentum have been 
earlier found for equation ([5]) in which only F^ u ^ or only ^ in papers [26J and 
[28] respectively. These solutions can be received from formula (127]) (the normalized 
stationary solutions from formula (130]) ). if either f 1 = or F^ v = . 



4. Discussion 



Let us discuss the physical meaning of the results obtained. For this purpose we shall 
consider vector and axial currents constructed with help of solution (1271) . The vector 
current is 

= V(x)j^(x) = (f/q°, (35) 

i.e. solution (127|) describes neutrino propagating with the velocity v = q/g°. At the 
same time the axial current is 

m 



A" = *(z)7Y*0*0 = Co'^. (36) 



Here 



where 



Q 

S» = -S? p (S S tp ) + [Sg + S? p (S S tp )] cos 29 - ±- e » vpX q v S Qp S tpX sin 29, (37) 



9 = (Nx)^f(^q) 2 - p 2 m 2 /m. (38) 
The spin vector £ can be expressed in terms of the four- vector components as 

< = S-n^-. (39) 
g u + m 

As the consequence of this fact we have 

C = C tp (CoC tp ) + [Co - Ct P (CoCt P )] cos 29 - [C tp x C ] sin 29. (40) 

Introduce a flight length L of a particle and an oscillations length L osc , using the 
relation 9 = nL/L osc . Since the scalar product (Nx) = r can be interpreted as the 
proper time of a particle, then the oscillation length is defined as 

Lose = (41) 

In this formula we use gaussian units and restore the neutrino magnetic moment /io- 

Hence if as a result of a certain process a neutrino arose with polarization £ , after 
travelling the distance L the probability for the neutrino to have polarization — £ is 
equal to 

W sf = [C xC tp ] 2 sm 2 (nL/L osc ). (42) 
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Consequently, if the condition (CoC p ) = is fulfilled, this probability can become unity, 
i.e. a resonance takes place. 

Thus, solution ( J27|) which is a linear combination of solutions ( 1301) describes a spin- 
coherent state of neutrino, in which a spin rotation takes place. Therefore, neutrino state 
with rotating spin is a pure state. Existence of such solutions is the direct consequence of 
the neutrino state description in terms of kinetic momentum. It should be stressed that 
as the result of calculations we obtained the complete system of neutrino wavef unctions, 
which show spin rotation properties. 

Such a spin behaviour has a simple quasi-classical interpretation. The 
antisymmetric tensor G^ v (see equation ( fTUl) ) can be written in the standard form 

G"" = (P,M), (43) 

where 

M = (/°q-g°f)/(2m), P = -[q x f]/(2m). (44) 

Vectors P and M are analogous to the polarization and the magnetization vectors of 
matter. It should be noted that the substitution F^ v =>- F^ u + G^ u implies that the 
magnetic H and electric D fields are shifted by the vectors M and P, respectively: 

H =^ B = H + M, D ^ E = D - P. (45) 

So restriction (jSJ) in the explicit form is 

(Ef) = 0, E/° - [B x f] = 0. (46) 

In the rest frame of the particle, equation (jSj) can be written down in the form 

C = 2[CxB ], (47) 

where the spin vector £ is related to four- vector by equation (1391) and the value B is 
the effective magnetic field in the neutrino rest frame which can be expressed in terms 
of quantities determined in the laboratory frame 

q(qB) 



B n = - 





1 


o — 






m 


1 








m 





q°B - [q x E] 



q° + m 



f 

2" 



[qxD]-^ + ^(7 -4^' n ' ^ 
[H J q° + m 2 V q° + m, 

Thus, the neutrino spin precesses around the direction B with the frequency u = 
2m|B |/g° = 27r|q|/(q ,0 L osc ). It is not difficult to see that the spin vector direction 
corresponding to stationary states Ct P is connected to the effective magnetic field as 
follows: 

B _ mL osc 

Ct " - - ^qT B °- (49) 

This fact explains the stationarity of states with Sq = S^ p . 

If neutrino possesses the fixed helicity in the initial state, i.e. 

^ = ^=i{|q| ) g°q/|q|}, Co = C sp = A (50) 
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then formula (I42p simplifies to the result widely discussed in the literature, which was 
mentioned in Introduction. 

Clearly, our arguments based on the consideration of vector and axial currents and 
the final formula (1421) are illustrations for neutrino propagation through background 
media. As it is commonly known [30J, in a relativistic case only the integrals of motion 
are well defined observables. 



5. Conclusions 



In this way we obtained the exact solutions of the Dirac-Pauli equation for neutrino 
in dense matter and electromagnetic field. It was demonstrated that if the neutrino 
production occurs in the presence of an external field and a dense matter, then its spin 
orientation is characterized by the vector S^ p instead of the vector S£ p . Using both 
the stationary ( 130]) and the nonstationary ( 1271) solutions it is possible to calculate the 
probabilities of various processes with neutrino in the framework of the Furry picture. 
Due to the time-energy uncertainty relation the considered states of neutrino can be 
generated only when the linear size of the area occupied by the electromagnetic field 
and the matter is comparable with the process formation length. This length is of the 
order of the oscillations length. 
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